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Scaling in long term data sets of geomagnetic indices and solar 
wind e as seen by WIND spacecraft. 
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We study scaling in fluctuations of the geomagnetic 
indices (AE, AU, and AL) that provide a measure of 
magnetospheric activity and of the e parameter which is 
a measure of the solar wind driver. Generalized structure 
function (GSF) analysis shows that fluctuations exhibit 
self-similar scaling up to about 1 hour for the AU in- 
dex and about 2 hours for AL, AE and e when the most 
extreme fluctuations over 10 standard deviations are ex- 
cluded. The scaling exponents of the GSF are found to be 
similar for the three AE indices, and to differ significantly 
from that of e. This is corroborated by direct comparison 
of their rescaled probability density functions. 



1. Introduction 

The statistical properties of fluctuations in geomag- 
netic indices and their relation to those in the solar 
wind, is a topic of considerable interest (see, e.g., [Sitnov 
et al, 2000; Tsurutani et al, 1990; Ukhorskiy et al, 2002; 
Voros et al., 1998]). Scaling has been identified as a key 
property of magnetospheric energy release in the form 
of bursty bulk flows in the magnetotail [Angelopoulos 
et al, 1992], "blobs" in the aurora [Lui et al., 2000], non- 
Gaussian fluctuations in geomagnetic indices [Hnat et al., 
2002, 2003a; Consolini et al., 1996] and in single station 
magnetometer data [Kovdcs et al., 2001; Voros et al, 
1998]. Models include Self-Organized Criticality (SOC) 
[Chang et al., 2003] (see also the review [Chapman and 
Watkins, 2001]) and multi-fractal models [Kovdcs et al, 
2001] related to those of turbulence[Consoimi et al., 1996; 
Voros et al, 1998]. 

These measures of scaling and non-Gaussian fluctua- 
tions in magnetospheric output need to be understood in 
the context of the system's driver, the solar wind, which is 
turbulent and thus also scaling. Other work has focussed 
on comparing properties of input parameters such as e 
and the indices (AE, AU and AL) to establish whether 
they are directly related. However, these studies have not 
provided a consistent answer. While Freeman et al. [2000] 
found that both the e and the AU and AL indices exhib- 
ited nearly identical scaling of burst lifetime probability 
density functions (PDFs), Uritsky et al. [2001] obtained 
quite different scalings for AE and the solar wind quan- 
tity v x B yz using spreading exponent methods motivated 
by SOC. Hnat et al. [2002, 2003a] used a PDF rescaling 
technique to characterize the fluctuation PDF of 4 years 
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e data from WIND and a 1 year data set of AE indices 
with fluctuations over a few standard deviations. Direct 
comparison of the PDF's functional form suggested close 
similarity to within statistical error. 

In this paper we use a larger 10-year data set for the 
AE indices to obtain a more accurate statistical determi- 
nation of the functional form of the PDF of fluctuations 
over a more extensive dynamic range, including charac- 
terization of extremal events up to 10 standard devia- 
tions for the first time. We apply structure functions to 
characterize and compare both the low and higher order 
moments for all quantities. A 4-year subset of the index 
data, corresponding to the same period in the solar cy- 
cle as that used to produce e, is used to facilitate this 
comparison. We then verify these results by direct ex- 
amination of the fluctuation PDF using the full 10-year 
AE indices dataset. 

2. Data Sets 

The AL, AU and AE index data sets investigated 
here comprise over 5.5 million, 1 minute averaged sam- 
ples from January 1978 to December 1988 inclusive. 
The e data set is identical to that used in Hnat et al. 
[2002, 2003a] and extends from January 1995 to Decem- 
ber 1998 inclusive. It includes intervals of slow and fast 
speed streams, e is defined (see [Hnat et al., 2002]) in 
SI units as e = «(B 2 //io)/o sin 4 (6/2), where lo ~ 7 Re 
and O = arct an ( | Bj, | /.Bz), and was calculated from the 
WIND spacecraft key parameter database [Lepping et al. , 
1995; Ogilvie et al, 1995]. The indices and e are from 
different time intervals and here we assume statistical 
stability over these long time intervals. 

3. Generalized Structure Functions 

Generalized structure functions (GSF), or generalized 
variograms, can be defined in terms of an average over 
time of a differenced variable Sx(t, r) = x(t + r) — x(t) as 
S m (T~) —< \8x(t, r)| m > [Rodriguez- Iturbe and Rinaldo, 
1997]. If Sx exhibits scaling with respect to r, then 
r C(™) a log-log plot of S m versus r should then reveal a 
straight line for each m with gradients ((m). If C( m ) = 
am (a constant) then the time series is self-similar with 
single scaling exponent a. 

In order to compare the scaling properties of the non- 
contemporaneous e and AE indices time series, we select 
a 4-year subinterval 1984 — 1987 from the AE indices at 
the same phase in the solar cycle as the e data. Fig- 
ure 1 shows the second order GSFs as measured by the 
standard deviations <t(t) = [S^t)] 1 / 2 of the fluctuation 
Sx(t,r). A scaling region is apparent between 2 7 and 
2 12 s where c(t) oc t h , where H is the Hurst exponent 
[£(2)/2]. The R goodness of fit analysis was performed 
to select the optimal power law region and gradient and 
results are summarized in Table 1. The upper limits of 
the scale regions T max are in good agreement with val- 
ues reported previously [Consolini and De Michelis, 1998; 
Takalo et al, 1993; Takalo and Timonen, 1998]. 

Any such single estimate of the H, whilst establish- 
ing the region of r over which there is scaling, does not 
fully characterize the properties of the time series. For 
example, a fractional Brownian motion (fBm) can be con- 
structed to share the same H value as AE, but the fBm 
series has Gaussian distributed increments Sx by defini- 
tion [Mandelbrot, 2002] whereas those of AE are non- 
Gaussian [Consolini and De Michelis, 1998; Hnat et al, 
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2002]. As discussed by Mandelbrot [2002] the similar val- 
ues arise because H aggregates two sources of scaling in 
monofractal random walks: persistence (the "Joseph" ef- 
fect) and heavy tails in the increments (the "Noah" ef- 
fect). In the above example the anomalous value of H 
for fBm comes just from the Joseph effect, whilst for AE 
the Noah effect must be at work. Furthermore, estimat- 
ing H by only one method may not distinguish a fractal 
time series from a discontinuous one [ Watkins et al, 2001; 
Katsev and L'Heureux, 2003]. We thus turn next to the 
higher order m values of (,{m). 

Figure 2 shows scaling exponents C,(m) derived from 
raw GSFs with m varying between —1 and 8 for the 8e 
and AE indices fluctuations. These suggest the depar- 
ture of higher orders from self-similarity, i.e., C(m) de- 
parts from a straight line. The inset of this figure shows 
the origin of these C( m ) values for 5AU and m = 1, 7. 
Only the first four orders exhibit clear linear behavior ex- 
pected in the scaling region. For higher orders, the value 
of £ very strongly depends on the assumed extent of the 
scaling region to which one fits a straight line. In princi- 
ple, C( m ) can be obtained for any m. However, errors do 
not contribute uniformly over m, for example, the largest 
fluctuations that affect large m, are statistically poorly 
resolved, whereas the smallest fluctuations (Sx — > 0) arc 
dominated by instrument thresholds. For the latter rea- 
son we will exclude m = — 1 for 8e as Se — > is not well 
determined through its definition. 

Conditioned GSFs quantify the impact of intermit- 
tency on fluctuations of different sizes by imposing a 
threshold A on the event size [Kovdcs et al., 2001]. Here, 
this threshold will be based on the standard deviation of 
the differenced time series for a given r, A(t) — 10ct(t). 
This procedure allows us to exclude rare extreme fluctua- 
tions with large statistical errors which, for large m, could 
lead to a spurious departure from self-similar behavior. 
Alternatively, conditioning with different thresholds es- 
timates a maximum size for the fluctuations for which 
self-similarity is still valid. 

Following conditioning, log- log plots of S^(t) show 
good correspondence with straight line fits, shown for 
SAU in the inset of figure 3. This power law dependence 
holds between times already obtained from the R 2 anal- 
ysis performed for <j(t). The main plot then shows C(m) 
obtained from the conditioned S^,(t). All lines in the 
figure were fitted for moments between —1 (0 for e) and 
6 and then extended to the entire range of data. Scaling 
exponents obtained from this technique were unchanged 
for thresholds A(r) between 6<r and 12cr. 

Firstly, our analysis suggests that the statistics of 
the fluctuations for all four quantities are self-similar for 
times between 2 and ~ 100 minutes and fluctuations of 
size 5x < 10<t(t). Secondly, the scaling exponent a in 
£(m) = am that characterize this self-similar behavior, 
are identical within errors for fluctuations in the AE in- 
dices but different to that in e at the la level. 

4. Rescaling of Fluctuation PDFs 

Scaling of the GSFs can be related to scaling properties 
of the fluctuation PDFs [Hnat et al, 2002, 2003a] using 
the generic, model-independent rescaling method (e.g. 
[Mantegna and Stanley, 1995; Hnat et al., 2003b]) based 
on the rescaling of the PDFs P(Sx, r) of 5x(t, r) on differ- 
ent time scales t. If a time series exhibits statistical self- 
similarity, a single argument representation of the PDF 
can be found that is given by P(Sx,t) — r~ a P s (5xT~ a ), 
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where a is the rescaling exponent. We now express S, 
using the fluctuations' PDF, P{8x, r) as follows: 



Expressing the integral in (1) in terms of rescaled vari- 
ables P s and 8x s = SxT~ a shows that the scaling ex- 
ponent C( m ) is a linear function of m, C( m ) = ma, for 
a statistically self-similar process, as suggested here by 
figure 3. 

The exponent a is ideally obtained from the scaling 
of the peaks of the PDF P(0,r). However, the finite ac- 
curacy of the measurement may discretize the amplitude 
leading to errors in the peak values. Table 1 gives all 
scaling exponents, obtained by different methods. These 
yield consistent values of a, to within the errors. We will 
use a from the scaling of a(r) versus t. If the fluctu- 
ations are statistically self-similar, as suggested by our 
GSF analysis, then the unsealed PDFs P(8x, r) should 
collapse onto a single curve P s (5x s ). We applied PDF 
rescaling to the fluctuation PDFs of all quantities and 
obtained satisfactory collapse of the curves within the 
scaling regions. The \ 2 test applied to all quantities re- 
vealed that, for the scaling regions given above, the col- 
lapsed curves lie within 5 — 7% error band. 

Figure 4 shows the re-scaled fluctuation PDFs for the 
indices alone for r ~ 15 min. The 8x variable has been 
normalized to the rescaled standard deviation a s (r ~ 
15mm.) of P s in each case to facilitate this comparison. 
The inset of this figure shows the comparison for AU, 
AE and — AL fluctuations and these PDFs are nearly 
identical. These results are consistent with conclusions 
of the GSF analysis at the la level. 

Figure 5 shows the normalized PDFs P a (Sx a ) for Sx — 
5e, SAE and r m 15 min overlaid on a single plot. We 
can clearly distinguish between the PDFs of the 5e and 
AE indices' fluctuations. We obtain the same result re- 
peating this comparison for several values of r, within the 
scaling range T ma x- We have also verified that the func- 
tional form of the PDF are insensitive to the solar cycle 
within errors. The use of a larger, 10 year data set for 
the indices has reduced statistical scatter and expanded 
the dynamic range of the considered fluctuations as com- 
pared to the analysis given in [Hnat et al., 2002, 2003a], 
and would lead us to draw the opposite conclusion, that 
on time scales less than ~ 1 hour the AE index amplitude 
fluctuations are not driven linearly by those of the solar 
wind. We would also conclude that the difference seen at 
the la level in the scaling of the e and the indices is sig- 
nificant, even though they agree at the 2a level [Freeman 
et al, 2000]. 

5. Summary 

In this paper we have addressed an open question of 
the possible connection between the scaling properties of 
fluctuations in the solar wind driver and those observed in 
global measures of magnetospheric dynamics. We applied 
two statistical methods, generalized structure functions 
and PDF rescaling, to study the scaling of fluctuations 
in the e parameter and the magnetospheric indices AU, 
AL and AE. We find that, statistically, fluctuations in all 
four quantities are approximately self-similar when their 
size is limited to ~ 10a. This self-similarity extends to 
~ 1 — 1.5 hours. The scaling exponents of the AE indices 
are close to each other and are appreciably different to 
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that of the e parameter. 

The fluctuation PDFs of the AE indices, unlike that 
of Se, are asymmetric. Direct comparison of the PDFs 
for the fluctuations in the AU , AE and — AL index in- 
dicates that they are nearly identical. Whilst the low 
frequency behavior of the solar wind and the indices may 
be well correlated [Tsurutani et al, 1990], here we have 
concluded that, on time scales smaller than 1 hour the 
properties of the fluctuations in the solar wind and the 
indices differ in both amplitude and persistence. If the 
underlying physical origin of the auroral scaling is turbu- 
lence, then different scaling behavior implies a different 
type of turbulence, i.e., different dimensionality /topology 
or different relevant physics [Frisch, 1995] . If the under- 
lying physics is SOC or similar [Chang et al, 2003] then 
similar conclusions would still be drawn (c.f. [Uritsky 
et al., 2001]). However, at this point we also can not 
rule out the possibility that the way in which the in- 
dices are constructed "burns" information still present 
in the magnetometer data about the solar wind scaling, 
here possibly by changing either or both of the degree of 
persistence (power spectral slope) and the heavy-tailed 
property (see [Edwards et al, 2001] for a related prelim- 
inary investigation). 
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Figure 1. Scaling of the standard deviation of the PDFs 
of: o-e, o-AU index, A-AL index and D-the AE index. 
The plots have been offset vertically for clarity. Error 
bars are estimated assuming Gaussian statistics for the 
binned data. 




Order m 

Figure 2. Dependence of the scaling exponent C(m) of 
the raw GSF on moment order m. Inset shows the GSF 
Sm versus time lag r for AU. 
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Figure 3. Dependence of the scaling exponent of the 
conditioned GSF on moment order. Inset shows the con- 
ditioned GSF versus time lag r for AU. 



Quantity 


a from P(0, r) 


a from <t(t) 


a from GSF 


Tmax [min\ 


e 




0.31 ±0.04 


0.25 ±0.04 


~ 100 


AE 


-0.47 ±0.03 


0.46 ±0.02 


0.40 ±0.02 


~ 100 


AU 


-0.46 ± 0.03 


0.41 ±0.02 


0.37 ±0.02 


~ 60 


AL 


-0.45 ± 0.03 


0.44 ± 0.02 


0.36 ±0.03 


~ 100 



Table 1. Scaling indices derived from P(0, r), <t(t) and GSF 
power laws. 
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8 x /a (t=1 5 min) 

Figure 4. Direct comparison between the fluctuation 
PDFs for AE{U), AU{o) and AL(A), again at r = 15 
minutes. Inset shows overlaid PDFs of AU, AE and — AL 
fluctuations. Error bars as in Figure 1. 
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Figure 5. Direct comparison, for the particular choice 
t — 15 minutes, of the fluctuation PDFs for e (o) and 
AE index (□). Error bars as in Figure 1. 



